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Introduction 

The last two decades have witnessed the recognition of a 
significant relationship between the autonomic nervous 
system and cardiovascular mortality, including sudden 
cardiac death 11 " 41 . Experimental evidence for an associ- 
ation between a propensity for lethal arrhythmias and 
signs of either increased sympathetic or reduced vagal 
activity has encouraged the development of quantitative 
markers of autonomic activity. 

Heart rate variability (HRV) represents one of 
the most promising such markers. The apparently easy 
derivation of this measure has popularized its use. As 
many commercial devices now provide automated 
measurement of HRV, the cardiologist has been pro- 
vided with a seemingly simple tool for both research and 
clinical studies 151 . However, the significance and meaning 
of the many different measures of HRV are more 
complex than generally appreciated and there is a 
potential for incorrect conclusions and for excessive or 
unfounded extrapolations. 

Recognition of these problems led the European 
Society of Cardiology and the North American Society 
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of Pacing and Electrophysiology to constitute a Task 
Force charged with the responsibility of developing 
appropriate standards. The specific goals of this Task 
Force were to: standardize nomenclature and develop 
definitions of terms; specify standard methods of 
measurement; define physiological and pathophysio- 
logical correlates; describe currently appropriate clinical 
applications, and identify areas for future research. 

In order to achieve these goals, the members of 
the Task Force were drawn from the fields of mathemat- 
ics, engineering, physiology, and clinical medicine. The 
standards and proposals offered in this text should 
not limit further development but. rather, should 
allow appropriate comparisons, promote circumspect 
interpretations, and lead to further progress in the field. 

The phenomenon that is the focus of this report 
is the oscillation in the interval between consecutive 
heart beats as well as the oscillations between consecu- 
tive instantaneous heart rates. 'Heart Rate Variability' 
has become the conventionally accepted term to describe 
variations of both instantaneous heart rate and RR 
intervals. In order to describe oscillation in consecutive 
cardiac cycles, other terms have been used in the litera- 
ture, for example cycle length variability, heart period 
variability, RR variabiuty and RR interval tachogram, 
and they more appropriately emphasize the fact that ft is 
the interval between consecutive beats that Is being 
analysed rather than the heart rate per se. However, 
these terms have not gained as wide acceptance as HRV, 
thus we will use the term HRV in this document. 



Background 

The clinical relevance of heart rate variability was first 
appreciated In 1965 when Hon and Lee 161 noted that fetal 
distress was preceded by alterations in ihterbeat intervals 
before any appreciable change occurred in the heart rate 
itself. Twenty years ago, Sayers and others focused 
attention on the existence of physiological rhythms 
imbedded in the beat-to-beat heart rate signal 17 ' 10 *. 
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During the 1970s, Ewing ef al [u] devised a number of 
simple bedside tests of short-term RR differences to 
detect autonomic neuropathy in diabetic patients. The 
association of higher risk of post-infarction mortality 
with reduced HRV was first shown by Wolf et a!, in 
1977 [i2} In Akselrod et al. introduced power 

spectral analysis of heart rate fluctuations to quantita- 
tively evaluate beat-to-beat cardiovascular control 11 31 . 

These frequency-domain analyses contributed to 
the understanding of the autonomic background of RR 
interval fluctuations In the heart rate record 114 - 151 . The 
clinical Importance of HRV became apparent in the late 
1980s when It was confirmed that HRV was a strong and 
independent predictor of mortality following an acute 
myocardial infarction 116 " 481 . With the availability of new, 
digital, high frequency, 24-h multi-channel electro- 
cardiographic recorders, HRV has the potential to 
provide additional valuable insight into physiological 
and pathological conditions and to enhance risk 
stratification. 



Measurement of heart rate variability 

Time domain methods 

Variations in heart rate may be evaluated by a number 
of methods. Perhaps the simplest to perform are the time 
domain measures. With these methods either the heart 
rate at any point in time or the intervals between 
successive normal complexes are determined. In a con- 
tinuous electrocardiographic (ECG) record, each QRS 
complex is detected, and the so-called normal-to-normal 
(NN) intervals (that is all intervals between adjacent 
QRS complexes resulting from sinus node depolariza- 
tions), or the instantaneous heart rate Is determined. 
Simple time-domain variables that can be calculated 
Include the mean NN interval, the mean heart rate, the 
difference between the longest and shortest NN interval, 
the difference between night and day heart rate, etc. 
Other time-domain measurements that can be used are 
variations in instantaneous heart rate secondary to 
respiration, tilt, Valsalva manoeuvre, or secondary to 
phenylephrine infusion. These differences can be de- 
scribed as either differences in heart rate or cycle length. 

Statistical methods 

From a series of instantaneous heart rates or cycle 
intervals, particularly those recorded over longer 
periods, traditionally 24 h, more complex statistical 
time-domain measures can be calculated. These may be 
divided into two classes, (a) those derived from direct 
measurements of the NN intervals or instantaneous 
heart rate, and (b) those derived from the differences 
between NN intervals. These variables may be derived 
from analysis of the total electrocardiographic recording 
or may be calculated using smaller segments of the 
recording period. The latter method allows comparison 
of HRV to be made during varying activities, e.g. rest, 
sleep, etc. 



The simplest variable to calculate is the standard 
deviation of the NN interval (SDNN). i.e. the square root 
of variance. Since variance is mathematically equal to 
total power of spectral analysis, SDNN reflects all the 
cyclic components responsible for variability in the 
period of recording. In many studies. SDNN is calcu- 
lated over a 24-h period and thus encompasses both 
short-term high frequency variations, as well as the 
lowest frequency components seen in a 24-h period. As 
the period of monitoring decreases, SDNN estimates 
shorter and shorter cycle lengths. It should also be noted 
that the total variance of HRV increases with the length 
of analysed recording 091 . Thus, on arbitrarily selected 
ECGs, SDNN is not a well defined statistical quantity 
because of its dependence on the length of recording 
period. Thus, in practice, it is inappropriate to compare 
SDNN measures obtained from recordings of different 
durations. However, durations of the recordings used to 
determine SDNN values (and similarly other HRV 
measures) should be standardized. As discussed further 
in this document, short-term 5-min recordings and 
nominal 24 h long-term recordings seem to be 
appropriate options. 

Other commonly used statistical variables calcu- 
lated from segments of the total monitoring period 
include SDANN, the standard deviation of the average 
NN interval calculated over short periods, usually 5 min, 
which is an estimate of the changes in heart rate due to 
cycles longer than 5 min, and the SDNN index, the mean 
of the 5-min standard deviation of the NN interval 
calculated over 24 h, which measures the variability due 
to cycles shorter than 5 min. 

The most commonly used measures derived from 
interval differences include RMSSD, the square root of 
the mean squared differences of successive NN intervals, 
NN50, the number of interval differences of successive 
NN intervals greater than 50 ms, and pNN50 the pro- 
portion derived by dividing NN50 by the total number 
of NN intervals. All these measurements of short-term 
variation estimate high frequency variations in heart 
rate and thus are highly correlated (Fig. 1). 

Geometrical methods 

The series of NN intervals can also be converted into a 
geometric pattern, such as the sample density distribu- 
tion of NN interval durations, sample density distribu- 
tion of differences between adjacent NN intervals, 
Lorenz plot of NN or RR intervals, etc., and a simple 
formula is used which judges the variability based on the 
geometric and/or graphic properties of the resulting 
pattern. Three general approaches are used in geometric 
methods: (a) a basic measurement of the geometric 
pattern (e.g. the width of the distribution histogram at 
the specified level) is converted into the measure of 
HRV, (b) the geometric pattern is interpolated by a 
mathematically defined shape (e.g. approximation of the 
distribution histogram by a triangle, or approximation 
of the differential histogram by an exponential curve) 
and then the parameters of this mathematical shape are 
used, and (c) the geometric shape is classified into several 



Eur Heart J, Vol. 17, March 1996 



356 Task Force 




100000 



10 000 - 



1000 - 



o 100- 




0.001 



0.1 1 
pNN50 (%) 



100 



Figure 1 Relationship between the RMSSD and pNN50 (a), and pNNSO and 
NN50 (b) measures of HRV assessed from 857 nominal 24-h Hotter tapes recorded 
in survivors of acute myocardial infarction prior to hospital discharge. The NN50 
measure used in panel (b) was normalized in respect to the length of the recording 
(Data of St. George's Post-infarction Research Survey Programme.) 



pattern-based categories which represent different 
classes of HRV (e.g. elliptic, linear and triangular shapes 
of Lorenz plots). Most geometric methods require the 
RR (or NN) interval sequence to be measured on or 
converted to a discrete scale which is not too fine or too 
coarse and which permits the construction of smoothed 
histograms. Most experience has been obtained with 
bins approximately 8 ms long (precisely 7-8125 ms- 
1/128 s) which corresponds to the precision of current 
commercial equipment. 

The HRV triangular index measurement is the 
integral of the density distribution (i.e. the number of all 



NN intervals) divided by the maximum of the density 
distribution. Using a measurement of NN intervals on a 
discrete scale, the measure is approximated by the value: 

(total number of NN intervals)/ 

(number of NN intervals in the modal bin) 

which is dependent on the length of the bin, i.e. on the 
precision of the discrete scale of measurement. Thus, if 
the discrete approximation of the measure is used with 
NN interval measurement on a scale different to the 
most frequent sampling of 128 Hz, the size of the bins 
should be quoted. The triangular interpolation of NN 
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Figure 2 To perform geometrical measures oa the NN interval 
histogram, the sample density distribution D is constructed which 
assigns the number of equally long NN intervals to each value of 
their lengths. The most frequent NN interval length Xk established, 
that is ?=D(A) is the maximum of the sample density distribution 
D. The HRV triangular index is the value obtained by dividing the 
area integral of D by the maximum Y. When constructing the 
distribution D with a discrete scale on the horizontal axis, the value 
is obtained according to die formula 

HRV index= (total number of all NN intervals)/ Y. 

For the computation of the TTNN measure, the values TV and M 
are established on the time axis and a multilinear function q 
constructed such that q(/)=0 for f< Wand Mand q(A)= Y, and 
such that the integral 



J;°°(D(t)-qW) 2 dt 



is the minimum among all selections of all values TV and M. The 
TTNN measure is expressed in ms and given by die formula 
TTNN=Af- N. 



interval histogram CTINN) is the baseline width of the 
distribution measured as a base of a triangle, approxi- 
mating the NN interval distribution (the minimum 
square difference is used to find such a triangle). Details 
of computing the HRV triangular index and TINN are 
shown in Fig. 2. Both these measures express overall 
HRV measured over 24 h and are more influenced by 
the lower than by the higher frequencies 117 *. Other 
geometric methods are still in the phase of exploration 
and explanation. 

The major advantage of geometric methods lies 
in their relative insensitivity to the analytical quality of 
the series of NN intervals 12 " 1 . The major disadvantage is 
the need for a reasonable number of NN intervals to 
construct the geometric pattern. In practice, recordings 
of at least 20 min (but preferably 24 h) should be used 
to ensure the correct performance of the geometric 
methods, i.e. the current geometric methods are in- 
appropriate to assess short-term changes in HRV. 

Summary and recommendations 

The variety of time-domain measures of HRV is sum- 
marized in Table 1 . Since many of the measures correlate 
closely with others, the following four are recommended 
for time-domain HRV assessment: SDNN (estimate of 
overall HRV): HRV triangular index (estimate of overall 
HRV); SDANN (estimate of long-term components of 



HRV). and RMSSD (estimate of short-term compo- 
nents of HRV). Two estimates of the overall HRV are 
recommended because the HRV triangular index 
permits only casual pre-processing of the ECG signal. 
The RMSSD method is preferred to pNN50 and NN50 
because it has better statistical properties. 

The methods expressing overall HRV and Its 
long- and short-term components cannot replace each 
other. The method selected should correspond to the 
aim of each study. Methods that might be recommended 
for clinical practices are summarized in the Section 
entitled Clinical use of heart rate variability. 

Distinction should be made between measures 
derived from direct measurements of NN intervals or 
instantaneous heart rate, and from the differences 
between NN intervals. 

It is inappropriate to compare time-domain 
measures, especially those expressing overall HRV, 
obtained from recordings of diflFerent durations. 

Other practical recommendations are listed in 
the Section on Recording requirements together with 
suggestions related to the frequency analysis of HRV. 

Frequency domain methods 

Various spectral methods 1231 for the analysis of the 
tachogram have been applied since the late 1960s. Power 
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Table 1 


Selected time-domain xnea 


wresofHRV 


Variable 


Units 


Description 






Statistical measures 



SDNN 


ms 


SDANN 


ms 


RMSSD 


ms 


SDNN Index 


ms 


SDSD 


ms 


NN50 count 




pNN50 


% 



Standard deviation of all NN intervals. 

Standard deviation of the averages of NN intervals in all 5 min segments of the entire recording. 
The square root of the mean of the sum of the squares of differences between adjacent NN 
intervals. 

Mean of the standard deviations of all NN intervals for all 5 min segments of the entire recording. 
Standard deviation of differences between adjacent NN intervals. 

Number of pairs of adjacent NN intervals differing by more than 50 ms in the entire recording. 
Three variants are possible counting all such NN intervals pairs or only pairs in which the first or 
the second interval is longer. 

NN50 count divided by the total number of all NN intervals. 



Geometric measures 



HRV triangular index Total number of all NN intervals divided by the height of the histogram of aH NN intervals 

-, measured on a discrete scale with bins of 7 8125 ms (1/128 s). (Details in Fig. 2) 

TINN ms Baseline width of the minimum square difference triangular interpolation of the highest peak of the 

histogram of all NN intervals (Details in Fig, 2.) 
Differential index ms Difference between the widths of the histogram of differences between adjacent NN Intervals 

measured at selected heights (eg. at the levels of 1000 and 10 000 samples) 12 ". 
Logarithmic index Coefficient p of the negative exponential curve k- which Is the best approximation of the 

histogram of absolute differences between adjacent NN intervals IZ2, . 



spectral density (PSD) analysis provides the basic infor- 
mation of how power (i.e. variance) distributes as a 
function of frequency. Independent of the method 
employed, only an estimate of the true PSD of the 
signals can be obtained by proper mathematical 
algorithms. 

Methods for the calculation of PSD may be 
generally classified as non-parametric and parametric. In 
most instances, both methods provide comparable 
results. The advantages of the non-parametric methods 
are; (a) the simplicity of the algorithm employed (Fast 
Fourier Transform — FFT — in most of the cases) and 
(b) the high processing speed, whilst the advantages of 
parametric methods are: (a) smoother spectral compo- 
nents which can be distinguished independently of pre- 
selected frequency bands, (b) easy post-processing of the 
spectrum with an automatic calculation of low and high 
frequency power components and easy identification of 
the central frequency of each component, and (c) an 
accurate estimation of PSD even on a small number of 
samples on which the signal is supposed to maintain 
stationarity. The basic disadvantage of parametric 
methods is the need to verify the suitability of the 
chosen model and its complexity (i.e. the order of the 
model). 

Spectra! components 

Short-term recordings Three main spectral components 
are distinguished in a spectrum calculated from short- 
term recordings of 2 to 5 min* 7,1013 ' 15 - 241 : very low fre- 
quency (VLF). low frequency (LF), and high frequency 
(HF) components. The distribution of the power and the 
central frequency of LF and HF are not fixed but may 
vary in relation to changes in autonomic modulations of 



the heart period 115 * 24 ' 251 . The physiological explanation 
of the VLF component is much less defined and the 
existence of a specific physiological process attributable 
to these heart period changes might even be questioned. 
The non-harmonic component which does not have 
coherent properties and which is affected by algorithms 
of baseline or trend removal is commonly accepted as a 
major constituent of VLF. Thus VLF assessed from 
short-term recordings (e.g. <5 min) is a dubious meas- 
ure and should be avoided when interpreting the PSD of 
short-term ECGs. 

Measurement of VLF, LF and HF power com- 
ponents is usually made in absolute values of power 
(ms 2 ), but LF and HF may also be measured In normal- 
ized units (n.u.) U5 MI which represent the relative value 
of each power component In proportion to the total 
power minus the VLF component. The representation of 
LF and HF in n.u. emphasizes the controlled and 
balanced behaviour of the two branches of the auto- 
nomic nervous system. Moreover, normalization tends 
to minimize the effect on the values of LF and HF 
components of the changes in total power (Fig. 3). 
Nevertheless. n.u. should always be quoted with abso- 
lute values of LF and HF power in order to describe in 
total the distribution of power in spectral components. 

Long-term recordings Spectral analysis may also be 
used to analyse the sequence of NN intervals in the 
entire 24-h period. The result then includes an ultra-low 
frequency component (ULF) t in addition to VLF, LF 
and HF components. The slope of the 24-h spectrum can 
also be assessed on a log-log scale by linear fitting the 
spectral values. Table 2 lists selected frequency-domain 
measures. 
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Figure 3 Spectral analysis (antoregressive model, order 
12) of RR interval variability in a healthy subject at rest 
and during 90 s head-up tilt At rest, two major components 
of similar power dee detectable at low and high frequen- 
cies. During tilt, the LF component becomes dominant 
but, as total variance is reduced, the absolute power of LF 
appears unchanged compared to rest Normalization pro- 
cedure leads to predominant LF and smaller HF compo- 
nents, which express the alteration of spectral components 
due to tilt The pie charts show the relative distribution 
together with the absolute power of the two components 
represented by the area. During rest, the total variance of 
the spectrum was 1201 ms 2 . and its VLF, LF, and HF 
components were 586 ms 2 , 310 ms 2 , and 302 ms 2 , respec- 
tively. Expressed in normalized units, the LF and HF were 
48-95 n.nu and 47-78 n.u., respectively. The LF/HF ratio 
was 1-02. During tilt, the total variance was 671 ms 2 , and 
its VLF, LF, and HF components were 265 ms 2 , 308 ms 2 , 
and 95 ms 2 , respectively. Expressed in normalized units, 
the LF and HF were 75*96 n.u. and 23-48 n. a., respec- 
tively. The LF/HF ratio was 3-34. Thus note that, for 
instance, the absolute power of the LF component was 
slightly decreased during tilt whilst the normalized units of 
LF were substantially increased 



The problem of stationarity' is frequently dis- 
cussed with long-term recordings. If mechanisms 
responsible for heart period modulations of a certain 
frequency remain unchanged during the whole period of 
recording, the corresponding frequency component of 
HRV may be used as a measure of these modulations. If 
the modulations are not stable, interpretation of the 
results of frequency analysis is less well defined. In 
particular, physiological mechanisms of heart period 
modulations responsible for LF and HF power compo- 
nents cannot be considered stationary during the 24-h 
period* 25 *. Thus, spectral analysis performed in the entire 
24-h period as well as spectral results obtained from 
shorter segments (e.g. 5 min) averaged over the entire 
24-h period (the LF and HF results of these two 
computations are not different p6 27) ) provide averages of 
the modulations attributable to the LF and HF compo- 
nents (Fig. 4). Such averages obscure detailed informa- 
tion about autonomic modulation of RR intervals 
available in shorter recordings 1251 . It should be remem- 
bered that the components of HRV provide measure- 
ments of the degree of autonomic modulations rather 




o 



Frequency (Hz) 

Figure 4 Example of an estimate of power spectral 
density obtained from the entire 24-h interval of a long- 
term Holler recording. Only the LF and HF components 
correspond to peaks of the spectrum while the VLF and 
ULF can be approximated by a line in mis plot with 
logarithmic scales on both axes. The slope of such a line 
is the a measure of HRV. 

than of the level of autonomic tone* 281 and averages of 
modulations do not represent an averaged level of tone. 

Technical requirements and recommendations 
Because of the important differences in the interpreta- 
tion of the results, the spectral analyses of short- and 
long-term electrocardiograms should always be strictly 
distinguished, as reported in Table 2. 

The analysed ECG signal should satisfy several 
requirements in order to obtain a reliable spectral esti- 
mation. Any departure from the following requirements 
may lead to unreproducible results that are difficult to 
interpret. 

In order to attribute individual spectral compo- 
nents to well defined physiological mechanisms, such 
mechanisms modulating the heart rate should not 
change during the recording. Transient physiological 
phenomena may perhaps be analysed by specific meth- 
ods which currently constitute a challenging research 
topic, but which are not yet ready to be used in applied 
research. To check the stability of the signal in terms of 
certain spectral components, traditional statistical tests 
may be employed 1291 . 

The sampling rate has to be properly chosen. A 
low sampling rate may produce a jitter in the estimation 
of the R wave fiducial point which alters the spectrum 
considerably. The optimal range is 250-500 Mz or per- 
haps even higher 1301 , while a lower sampling rate (in any 
case ^100 Hz) may behave satisfactorily only if an 
algorithm of interpolation (e.g. parabolic) is used to 
refine the R wave fiducial point 01 - 321 . 

Baseline and trend removal (if used) may affect 
the lower components in the spectrum. It is advisable to 
check the frequency response of the filter or the behav- 
iour of the regression algorithm and to verify that the 
spectra] components of interest are not significantly 
affected. 
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Table 2 Selected frequency domain measures ofHRV 



Variable 


Units 


Description 
Analysis of short-term recordings (5 min) 


Frequency range 


5 min total power 


ms z 


The variance of NN intervals over the 


approximately £0-4 Hz 




ms 2 


temporal scg^ou 


VLF 


Power in very low frequency range 


£0.04 Hz 


LF 


ms 2 


Power in low frequency range 


0 04-015 Hz 


LF norm 


au. 


LF power in normalised units 








LF/(Total Power-VLF) x 100 




HF 


ins 2 


Power in high frequency range 


0-15-0-4 Hz 


HF norm 


n.iL 


HF power in normalised units 








HF/ (Total Power-VLF) x 100 




LF/HF 




Ratio LF [ms^VHF Ims 2 ] 








Analysis of entire 24 h 





Total power, ras 2 Variance of all NN intervals 

ULF ; \ ms 2 Power in the ultra low frequency range 

VLF ms 2 Power in the very low frequency range 

LF ms 2 Power in the low frequency range 

HF ms 2 Power In the high frequency range 

a Slope of the linear interpolation of the 
spectrum In a log-log scale 



approximately <0-4 Hz 
£0-003 Hz 
0-003-0-04 Hz 
004-015 Hz 
015-0-4 Hz 

approximately £0-04 Hz 



The choice of QRS fiducial point may be critical. 
It is necessary to use a well tested algorithm (i.e. 
derivative+threshold, template, correlation method, 
etc) in order to locate a stable and noise-independent 
reference point 1331 . A fiducial point localized far within 
the QRS complex may also be influenced by varying 
ventricular conduction disturbances. 

Ectopic beats, arrhythmic events, missing data 
and noise effects may alter the estimation of the PSD of 
HRV. Proper interpolation (or linear regression or simi- 
lar algorithms) on preceding/successive beats on the 
HRV signals or on its autocorrelation function may 
reduce this error. Preferentially, short-term recordings 
which are free of ectopy, missing data, and noise should 
be used. In some circumstances, however, acceptance of 
only ectopic-free short-term recordings may introduce 
significant selection bias. In such cases, proper inter- 
polation should be used and the possibility of the results 
being influenced by ectopy should be considered 1341 . The 
relative number and relative duration of RR intervals 
which were omitted and interpolated should also be 
quoted. 

Algorithmic standards and recommendations 
The series of data subjected to spectral analysis can be 
obtained in different ways. A useful pictorial represen- 
tation of the data Is the discrete event series (DES), that 
is the plot of Rj-Ri-j interval vs time (indicated at Ri 
occurrence) which is an irregularly time-sampled signal. 
Nevertheless, spectral analysis of the sequence of 
instantaneous heart rates has also been used in many 
studies 1261 . 

The spectrum of the HRV signal is generally 
calculated either from the RR interval tachogram (RR 
durations vs number of progressive beats — see Fig. 



5a.b) or by interpolating the DES, thus obtaining a 
continuous signal as a function of time, or by calculating 
the spectrum of the counts — unitary pulses as a func- 
tion of time corresponding to each recognised QRS 
complex* 351 . Such a choice may have implications on the 
morphology, the measurement units of the spectra and 
the measurement of the relevant spectral parameters. In 
order to standardize the methods, the use of the RR 
interval tachogram with the parametric method, or the 
regularly sampled interpolation of DES with the non- 
parametric method may be suggested; nevertheless, 
regularly sampled interpolation of DES is also suitable 
for parametric methods. The sampling frequency of 
interpolation of DES has to be sufficiently high that the 
Nyquist frequency of the spectrum is not within the 
frequency range of interest 

Standards for non-parametric methods (based 
upon the FFT algorithm) should include the values 
reported in Table 2. the formula of DES interpolation, 
the frequency of sampling the DES interpolation, the 
number of samples used for the spectrum calculation, 
and the spectral window employed (Hann, Harriming, 
and triangular windows are most frequently used) 13 ". 
The method of calculating the power in respect of the 
window should also be quoted. In addition to require- 
ments described in other parts of this document, each 
study employing the non-parametric spectral analysis of 
HRV should quote all these parameters. 

Standards for parametric methods should 
Include the values reported in Table 2, the type of the 
model used, the number of samples* the central fre- 
quency for each spectral component (LF and HF) and 
the value of the model order (numbers of parameters). 
Furthermore, statistical figures have to be calculated in 
order to test the reliability of the model. The prediction 



Eur Heart J, Vol. 17, March 1996 



Standards of heart rate variability 361 



Rest 
tachogram 



1.0 



0.5 
0.4 







7 ii = 256 




7 Mean - 842.5 ms 
i 1 


O 2 = 1784 ms 2 

I 



100 200 
Beat* 



1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 



Tilt 
tachogram 



"(b) 




. n = 256 




[ Mean ~ 564.7 ms 


a 2 » 723 ms 2 


. ! 





100 200 
Beat* 



0.015 



0.010 - 



I 

g 



0.005 



(c) 





Frequency 


Power 
ms 2 


Power 
n.u. 


VLF- 


0.00 


785 




LF ; 


0.11 


479 


47.95 


HF 


O.Z4 


450 


45.05 



VLF 



LF/HF = 1.06 
PEWT>3 
OOT = 9 
p = 9 




0.015 



0.010 
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(cO 





Frequency 


Power 
ms* 


Power 
n.u. 


VLF 


0.00 


192 




LF 


0.09 


413 


77,78 


HF 


0.24 


107 


20.15 



0.1 0.2 0.3 0.4 
Frequency (Hz) 




LF/HF = 3.6 
PEWT> 11 
OOT= 15 
p = 15 
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Figures Interval tachogram of 256 consecutive RR values in a normal subject at supine rest (a) and after 
head-up t ill (b ). The HRV spectra are shown, calculated by parametric autoregressive modelling (c and d), 
and by a FFT based non-parametric algorithm (e and f). Mean values (m), variances (s 2 ) and the* number 
(n) of samples are indicated. For (c) and (d), VLF, LF and HF central frequency, power in absolute value 
and power in normalized units (n.u,) are also indicated together with the order p of the chosen model and 
minimal values of PEYVT and OOT which satisfy the tests. In (e) and (f) , the peak frequency and the power 
of VLF, LF and HF were calculated by integrating the PSD in die defined frequency bands. The window 
type is also specified. In panels (c) to (f), the LF component is indicated by dark shaded areas and the HF 
component by light shaded areas. 
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Table 3 Approximate correspondence of time domain 
and frequency domain methods applied to 24-b ECG 



Time domain variable 


Approximate frequency 
domain correlate 


SDNN 


Total power 


HRV triangular index 


Total power 


TINN 


Total power 


SDANN 


ULF 


SDNN index 


Mean of 5 rain total power 


RMSSD 


HF 


SDSD 


HF 


NN50 count 


HF 


pNN50 


HF 


Differential index 


HF 


Logarithmic index 


HF 



error whiteness test (PEWT) provides information about 
the 'goodness' of the fitting model* 371 while the optimal 
order test (OOI) checks the suitability of the order of 
the model used 1381 . There are different possibilities of 
performing OOT which include final prediction error 
and Akaike information criteria. The following opera- 
tive criterion for choosing the order p of an autoregres- 
sive model might be proposed: the order shall be in the 
range 8-20. fulfilling the PEWT test and complying with 
the OOT test (p w min(OOT)) . 

Correlation and differences between time and frequency 
domain measures 

When analysing stationary short-term recordings, more 
experience and theoretical knowledge exists on the 
physiological interpretation of the frequency-domain 
measures compared to the time-domain measures 
derived from the same recordings. 



However, many time- and frequency-domain 
variables measured over the entire 24-h period are 
strongly correlated with each other (Table 3). These 
strong correlations exist because of both mathematical 
and physiological relationships. In addition, the physio- 
logical interpretation of the spectral components calcu- 
lated over 24 h is difficult, for the reasons mentioned 
(section entitled Long-term recordings). Thus, unless 
special investigations are performed which use the 24-h 
HRV signal to extract information other than the usual 
frequency components (e.g. the log-log slope of spectro- 
gram), the results of frequency-domain analysis are 
equivalent to those of time-domain analysis, which is 
easier to perform. 



Rhythm pattern analysis 

As illustrated in Fig. G m , the time-domain and spectral 
methods share limitations imposed by the irregularity of 
the RR series. Clearly different profiles analysed by these 
techniques may give identical results. Trends of decreas- 
ing or increasing cycle length are in reality not symmetri- 
cal 140 - 411 as heart rate accelerations are usually followed 
by a faster decrease. In spectral results, this tends to 
reduce the peak at the fundamental frequency, and to 
enlarge its basis. This leads to the idea of measuring 
blocks of RR intervals detennined by properties of the 
rhythm and investigating the relationship of such blocks 
without considering the internal variability. 

Approaches derived from the time-domain and 
the frequency-domain have been proposed in order to 
reduce these difficulties. The interval spectrum and spec- 
trum of counts methods lead to equivalent results (d. 
Fig. 6) and are well suited to investigate the relationship 
between HRV and the variability of other physiological 
measures. The interval spectrum is well adapted to link 
RR intervals to variables defined on a beat-to-beat basis 



(a) 



(b) 



(d) 



500 



1000 



Figure 6 Example of four synthesised time series with identical means, 
standard deviations, and ranges. Series (c) and (d) also have identical 
autocorrelation functions and therefore identical power spectra. Reprinted 
with permission 1391 . 
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(e.g. blood pressure). The spectrum of counts is prefer- 
able if RR intervals are related to a continuous signal 
(e.g. respiration), or to the occurrence of special events 
(e.g. arrhythmia). 

The "peak-valley* procedures are based either on 
the detection of the summit and the nadir of oscilla- 
tions 142,431 or on the detection of trends of heart rate 1441 . 
The detection may be limited to short-term changes 1421 
but it can be extended to longer variations: second and 
third order peaks and troughs* 431 or stepwise increase of 
a sequence of consecutive increasing or decreasing cycles 
surrounded by opposite trends 1441 . The various oscilla- 
tions can be characterized on the basis of the heart 
rate accelerating or slowing, the wavelength and/or the 
amplitude. In a majority of short- to mid-term record- 
ings, the results are correlated with frequency compo- 
nents of HRV I4 ^wThe correlations, however, tend to 
diminish as the wavelength of the oscillations and the 
recording duration increase. Complex demodulation 
uses the techniques of interpolation and detr ending 14 61 
and provides the time resolution necessary to detect 
short-term heart rate changes, as well as to describe 
the amplitude and phase of particular frequency 
components as functions of time. 



Non-linear methods 

Non-linear phenomena are certainly involved in the 
genesis of HRV. They are determined by complex inter- 
actions of haemodynamic. electrophysiological and 
humoral variables, as well as by autonomic and central 
nervous regulations. It has been speculated that analysis 
of HRV based on the methods of non-linear dynamics 
might elicit valuable information for the physiological 
interpretation of HRV and for the assessment of the risk 
of sudden death. The parameters which have been used 
to measure non-linear properties of HRV include 1/f 
scaling of Fourier spectra 147191 , H scaling exponent, and 
Coarse Graining Spectral Analysis (CGSA)* 48 *. For data 
representation, Poincare sections, low-dimension attrac- 
ted plots, singular value decomposition, and attractor 
trajectories have been used. For other quantitative 
descriptions, the D 2 correlation dimension, Lyapunov 
exponents, and Kolmogorov entropy have been 
employed 1495 . 

Although in principle these techniques have been 
shown to be powerful tools for characterization of 
various complex systems, no major breakthrough has 
yet been achieved by their application to bio-medical 
data including HRV analysis. It is possible that integral 
complexity measures are not adequate to analyse bio- 
logical systems and thus, are too insensitive to detect the 
non-linear perturbations of RR interval which would be 
of physiological or practical importance. More encour- 
aging results have been obtained using differential, 
rather than Integral complexity measures, e.g. the scaling 
index method 1 - 511 . However, no systematic study has 
been conducted to investigate large patient populations 
using these methods. 



At present, the non-linear methods represent 
potentially promising . tools for HRV assessment, but 
standards are lacking and the full scope of these meth- 
ods cannot be assessed. Advances in technology and the 
interpretation of the results of non-linear methods are 
needed before these methods are ready for physiological 
and clinical studies. 

Stability and reproducibility of HRV 
measurement 

Multiple studies have demonstrated that short-term 
measures of HRV rapidly return to baseline after tran- 
sient perturbations Induced by such manipulations as 
mild exercise, administration of short acting vasodila- 
tors, transient coronary occlusion, etc. More powerful 
stimuli, such as maximum exercise or administration of 
long acting drugs may result in a much more prolonged 
interval before return to control values. 

There are far fewer data on the stability of 
long-term measures of HRV obtained from 24-h ambu- 
latory monitoring. Nonetheless, the same amount of 
data available suggest great stability of HRV measures 
derived from 24-h ambulatory monitoring in both 
normal subjects 152 * 531 and in the post-infarction 1541 and 
ventricular arrhythmia 1551 populations. There also exists 
some fragmentary data to suggest that stability of HRV 
measures may persist for months and years. Because 
24-h Indices seem to be stable and free of placebo 
effect, they may be ideal variables with which to assess 
intervention therapies. 



Recording requirements 

ECG signal 

The fiducial point recognised on the ECG tracing which 
identifies a QRS complex may be based on the maxi- 
mum or baricentrum of the complex, on the determina- 
tion of the maximum of an interpolating curve, or found 
by matching with a template or other event markers. In 
order to localize the fiducial point, voluntary standards 
for diagnostic ECG equipment are satisfactory in 
terms of signal/noise ratio, common mode rejection, 
bandwidth, etc.* 56 ' An upper-band frequency cut-off 
substantially lower than that established for diagnostic 
equipment (-200 Hz) may create a jitter in the recogni- 
tion of the QRS complex fiducial point, introducing 
an error of measured RR intervals. Similarly, limited 
sampling rate induces an error in the HRV spectrum 
which increases with frequency, thus affecting more high 
frequency components 1311 . An interpolation of the 
undersampled ECG signal may decrease this error. With 
proper interpolation, even a 100 Hz sampling rate can be 
sufficient 1321 . 

When using solid-state storage recorders, data 
compression techniques have to be carefully considered 
in terms of both the effective sampling rate and the 



Eur Heart J. Vol. 17. March 1996 



364 Task Force 



quality of reconstruction methods which may yield 
amplitude and phase distortion 1571 . 

Duration and circumstances ofECG recording 
In studies researching HRV, the duration of recording is 
dictated by the nature of each investigation. Standardiz- 
ation is needed, particularly in studies investigating the 
physiological and clinical potential of HRV. 

Frequency-domain methods should be preferred 
to the time-domain methods when investigating short- 
term recordings. The recording should last for at least 10 
times the wavelength of the lower frequency bound of 
the investigated component, and, in order to ensure 
the stability of the signal, should not be substantially 
extended. Thus, * recording of approximately 1 min is 
needed to assess the HF components of HRV while 
approximately 2 min are needed to address the LF 
component. In order to standardize different studies 
investigating short-term HRV, 5 min recordings of a 
stationary system are preferred unless the nature of the 
study dictates another design. 

Averaging of spectral components obtained from 
sequential periods of time is able to minimize the error 
imposed by the analysis of very short segments. Never- 
theless, if the nature and degree of physiological heart 
period modulations changes from one short segment of 
the recording to another, the physiological interpreta- 
tion of such averaged spectral components suffers from 
the same intrinsic problems as that of the spectral 
analysis of long-term recordings and warrants further 
elucidation. A display of stacked series of sequential 
power spectra (e.g. over 20 min) may help confirm 
steady state conditions for a given physiological state. 

Although the time-domain methods, especially 
the SDNN and RMSSD methods, can be used to 
investigate recordings of short durations, the frequency 
methods are usually able to provide more easily inter- 
pretable results in terms of physiological regulations. In 
general, the time-domain methods are Ideal for the 
analysis of long-term recordings (the lower stability of 
heart rate modulations during long-term recordings 
makes the results of frequency methods less easily inter- 
pretable). The experience shows that a substantial part 
of the long-term HRV value is contributed by the 
day-night differences. Thus the long-term recording 
analysed by the time-domain methods should contain at 
least 18 h of analysable ECG data that includes the 
whole night. 

Little is known about the effects of the environ- 
ment (e.g. type and nature of physical activity and 
of emotional circumstances) during long-term ECG 
recordings. For some experimental designs, environ- 
mental variables should be controlled and in each study, 
the character of the environment should always be 
described. The design of investigations should also 
ensure that the recording environment of individual 
subjects is similar. In physiological studies comparing 
HRV in different well-defined groups, the differences 
between underlying heart rate should also be properly 
acknowledged. 



Editing of the RR interval sequence 
The errors imposed by the imprecision of the NN 
interval sequence are known to affect substantially the 
results of statistical time-domain and all frequency- 
domain methods. It is known that casual editing of the 
RR interval data is sufficient for the approximate assess- 
ment of total HRV by the geometric methods, but it is 
not known how precise the editing should be to ensure 
correct results from other methods. Thus when using the 
statistical time-domain and/or frequency-domain meth- 
ods, the manual editing of the RR data should be 
performed to a very high standard ensuring correct 
identification and classification of every QRS complex. 
Automatic 'filters' which exclude, some intervals from 
the original RR sequence (e.g. those differing by more 
than 20% from the previous interval) should not replace 
manual editing as they are known to behave unsatisfac- 
torily and to have undesirable effects leading potentially 
to errors 1581 . 

Suggestions for standardisation of commercial 
equipment 

Standard measurement of HRV Commercial equipment 
designed to analyse short-term HRV should incorporate 
non-parametric and preferably also parametric spectra] 
analysis. In order to iruniinize the possible confusion 
imposed by reporting the components of the cardiac 
beat-based analysis in time-frequency components, the 
analysis based on regular sampling of the tachograms 
should be offered in all cases. Non-parametric spectral 
analysis should employ at least 512 but preferably 1024 
points for 5 min recordings. 

Equipment designed to analyse HRV in long- 
term recordings should implement time-domain 
methods including all four standard measures (SDNN, 
SDANN, RMSSD, and HRV triangular index). In 
addition to other options, the frequency analysis should 
be performed in 5 min segments (using the same preci- 
sion as with the analysis of short-term ECGs). When 
performing the spectral analysis of the total nominal 
24-h record in order to compute the whole range of HF, 
LF, VLF and ULF components, the analysis should 
be performed with a similar precision of periodogram 
sampling, as suggested for the short-term analysis, e.g. 
using 2 18 points. 

The strategy of obtaining the data for the HRV 
analysis should copy the design outlined in Fig. 7. 

Precision and testing of commercial equipment In order 
to ensure the quality of different equipment involved 
in HRV analysis and to find an appropriate balance 
between the precision essential to research and clinical 
studies and the cost of the equipment required, indepen- 
dent testing of all equipment is needed. As the potential 
errors of the HRV assessment include inaccuracies in the 
identification of fiducial points of QRS complexes, the 
testing should include all the recording, replay, and 
analysis phases. Thus, it seems ideal to test various 
equipment with signals (e.g. computer simulated) of 
known HRV properties rather than with existing 
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Figure 7 How chart summarizing individual steps used 
when recording and processing the ECG signal in order to 
obtain data for H$V analysis. 



databases of already digitized ECGs. When employing 
commercial equipment in studies investigating physi- 
ological and clinical aspects of HRV, independent tests 
of the equipment used should always be required. A 
possible strategy for testing of commercial equipment is 
proposed in Appendix B. Voluntary industrial standards 
should be developed adopting this or similar strategy. 



In order to minimize the errors caused by improperly 
designed or incorrectly used techniques, the following 
points are recommended: 

The ECG equipment used should satisfy the 
current voluntary industrial standards in terms of signal/ 
noise ratio, common mode rejection, bandwidth, etc. 

Solid-state recorders used should allow signal 
reconstruction without amplitude and phase distortion; 
long-term ECG recorders using analogue magnetic 
media should accompany the signal with phase-locked 
time tracking. 

Commercial equipment used to assess HRV 
should satisfy the technical requirements listed in the 
Section on Standard Measurement of HRV and its 
performance should be independently tested. 

In order to standardize physiological and clinical 
studies, two types of recordings should be used whenever 
possible: (a) short-term recordings of 5 mln made 
under physiologically stable conditions processed by 
fr equency-domain methods, and/or (b) nominal 24-h 
recordings processed by time-domain methods. 

When using long-term ECGs in clinical studies, 
individual, subjects should be recorded under fairly 
similar conditions and in a fairly similar environment. 

When using statistical time--domain or 
frequency-domain methods, the complete signal should 
be carefully edited using visual checks and manual 
corrections of individual RR intervals and QRS complex 
classifications. Automatic 'filters' based on hypotheses 
on the logic of RR interval sequence (e.g. exclusion of 



RR intervals according to a certain prematurity thresh- 
old) should not be relied on when ensuring the quality of 
the RR interval sequence. 



Physiological correlates of heart rate 
variability 

Physiological correlates of HRV components 

Autonomic influences of heart rate 
Although cardiac automaticity is intrinsic to various 
pacemaker tissues, heart rate and rhythm are largely 
under the control of the autonomic nervous system* 59 *. 
The parasympathetic Influence on heart rate is mediated 
via release of acetylcholine by the vagus nerve. 
Muscarinic acetylcholine receptors respond to this 
release mostly by an increase in cell membrane K+ 
conductance 16 *" 62 *. Acetylcholine also inhibits the 
hyperpolarizatlon-activated "pacemaker* current If* 63 * 641 . 
The Tk decay' hypothesis 165 * proposes that pacemaker 
depolarization results from slow deactivation of the 
delayed rectifier current, Ik, which, due to a time- 
independent background inward current, causes diasto- 
lic depolarization* 65 - 661 . Conversely, the 'If activation* 
hypothesis 1671 suggest that following action potential 
termination. If provides a slowly activating inward cur- 
rent predominating over decaying Ik, thus initiating 
slow diastolic depolarization. 

The sympathetic influence on heart rate is medi- 
ated by release of epinephrine and norepinephrine. 
Activation of ^-adrenergic receptors results In cyclic 
AMP mediated phosphorylation of membrane proteins 
and increases in ICaJJ 681 and in If 169 * 701 . The end result is 
an acceleration of the slow diastolic depolarization. 

Under resting conditions, vagal tone prevails 171 * 1 
and variations in heart period are largely dependent on 
vagal modulation 1721 . The vagal and sympathetic activity 
constantly interact. As the sinus node is rich in acetyl- 
cholinesterase, the effect of any vagal impulse is brief 
because the acetylcholine is rapidly hydrorysed. Para- 
sympathetic influences exceed sympathetic effects 
probably via two independent mechanisms: a cholin- 
ergicaUy induced reduction of norepinephrine released 
in response to sympathetic activity, and a cholinergic 
attenuation of the response to a adrenergic stimulus. 

Components of HRV 

The RR interval variations present during resting con- 
ditions represent a fine tuning of the beat-to-beat control 
mechanisms^ 3 * 713 . Vagal afferent stimulation leads to 
reflex excitation of vagal efferent activity and inhibition 
of sympathetic efferent activity 175 *. The opposite reflex 
effects are mediated by the stimulation of sympathetic 
afferent activity 1761 . Efferent vagal activity also appears 
to be under 'tonic* restraint by cardiac afferent sympa- 
thetic activity 071 . Efferent sympathetic and vagal activi- 
ties directed to the sinus node are characterized by 
discharge largely synchronous with each cardiac cycle 
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which can be modulated by central (e.g. vasomotor and 
respiratory centres) and peripheral (e.g. oscillation in 
arterial pressure and respiratory movements) oscilla- 
tors {241 . These oscillators generate rhythmic fluctuations 
in efferent neural discharge which manifest as short and 
long-term oscillation in the heart period. Analysis of 
these rhythms may permit inferences on the state and 
function of (a) the central oscillators, (b) the sympa- 
thetic and vagal efferent activity, (c) humoral factors, 
and (d) the sinus node. 

An understanding of the modulatory effects of 
neural mechanisms on the sinus node has been enhanced 
by spectral analysis of HRV. The efferent vagal activity 
is a major contributor to the HF component, as seen in 
clinical and experimental observations of autonomic 
manoeuvres such as electrical vagal stimulation, mus- 
carinic receptor blockade, and vagotomy 13 ' 14,24 *. More 
controversial is the interpretation of the LF component 
which is considered by some 124 * 78-801 as a marker of 
sympathetic modulation (especially when expressing it in 
normalized units) and by others 113,811 as a parameter that 
includes both sympathetic and vagal influences. This 
discrepancy is due to the fact that in some conditions, 
associated with sympathetic excitation, a decrease in the 
absolute power of the LF component is observed. It is 
important to recall that during sympathetic activation 
the resulting tachycardia is usually accompanied by a 
marked reduction in total power, whereas the reverse 
occurs during vagal activation. When the spectral com- 
ponents are expressed in absolute units (ms 2 ). the 
changes in total power influence LF and HF in the same 
direction and prevent the appreciation of the fractional 
distribution of the energy. This explains why in supine 
subjects under controlled respiration atropine reduces 
both LF and HF 1141 and why during exercise LF is 
markedly reduced 1241 . This concept is exemplified in 
Fig. 3 showing the spectral analysis of HRV in a normal 
subject during control supine conditions and 90* 
head-up tilt. Due to the reduction in total power, 
LF appears as unchanged if considered in absolute 
units. However, after normalization an increase in LF 
becomes evident. Similar results apply to the LF/HF 
ratio 1821 . 

Spectral analysis of 24-h recordings 124 " 251 shows 
that in normal subjects LF and HF expressed in normal- 
ized units exhibit a circadian pattern and reciprocal 
fluctuations, with higher values of LF in the daytime and 
of HF at night. These patterns become undetectable 
when a single spectrum of the entire 24-h period is used 
or when spectra of subsequent shorter segments are 
averaged. In long-term recordings, the HF and LF 
components account for approximately 5% of total 
power. Although the ULF and VLF components 
account for the remaining 95% of total power, their 
physiological correlates are still unknown. 

LF and HF can increase under different condi- 
tions. An increased LF (expressed in normalized units) is 
observed during 90* tilt, standing, mental stress and 
moderate exercise in healthy subjects, and during mod- 
erate hypotension, physical activity and occlusion of a 



coronary artery or common carotid arteries in conscious 
dogs* 24 ™. Conversely, an increase in HF is induced by 
controlled respiration, cold stimulation of the face and 
rotational staidi 124 - 781 . 

Summary and recommendations for interpretation of 
HRV components 

Vagal activity is the major contributor to the HF 
component. 

Disagreement exists in respect of the LF compo- 
nent. Some studies suggest that LF, when expressed in 
normalized units, is a quantitative marker for sympa- 
thetic modulations, other studies view LF as reflecting 
both sympathetic and vagal activity. Consequendy, the 
LF/HF ratio is considered by some investigators to 
mirror sympatho/vagal balance or to reflect sympathetic 
modulations. 

Physiological interpretation of lower frequency 
components of HRV (that is of the VLF and ULF 
components) warrants further elucidation. 

It is important to note that HRV measures 
fluctuations in autonomic inputs to the heart rather than 
the mean level of autonomic inputs. Thus both auto- 
nomic withdrawal and a saturatingly high level of 
sympathetic input leads to diminished HRV* 281 . 

Changes of HRV related to specific 
pathologies 

A reduction of HRV has been reported in several 
cardiological and non -cardiological diseases 124,78 ' 81,831 . 

Myocardial infarction 

Depressed HRV after MI may reflect a decrease in vagal 
activity directed to the heart which leads to prevalence 
of sympathetic mechanisms and to cardiac electrical 
instability. In the acute phase of MI, the reduction in 
24-h SDNN is significantly related to left ventricular 
dysfunction, peak creatine kinase, and Killip class 1841 . 

The mechanism by which HRV is transiently 
reduced after MI and by which a depressed HRV is 
predictive of the neural response to acute MI is not yet 
defined, but it is likely to involve derangements in the 
neural activity of cardiac origin. One hypothesis 1851 
Involves cardio-cardiac sympatho-syinpathetic 188 ' 871 and 
sympatho-vagal reflexes^ 1 and suggests that the changes 
in the geometry of a beating heart due to necrotic and 
non-contracting segments may abnormally increase the 
firing of sympathetic afferent fibres by mechanical dis- 
tortion of the sensory ending ^7B,87<88, . This sympathetic 
excitation attenuates the activity of vagal fibres directed 
to the sinus node. Another explanation, especially 
applicable to marked reduction of HRV, is the 
reduced responsiveness of sinus nodal cells to neural 
modulations 182 - 851 . 

Spectral analysis of HRV in patients surviving an 
acute MI revealed a reduction in total and in the 
individual power of spectral components 1891 . However, 



Eur Heart J. Vol. 17. March 1996 



Standards of heart rate variability 367 



when the power of LF and HF was calculated in 
normalized units, an increased LF and a diminished HF 
were observed during both resting controlled conditions 
and 24-h recordings analysed over multiple 5 mln 
periods 190311 . These changes may indicate a shift of 
sympatho-vagal balance towards sympathetic predomi- 
nance and reduced vagal tone. Similar conclusions were 
obtained by considering the changes in the LF/HF ratio. 
Hie presence of an alteration in neural control mecha- 
nisms was also reflected by the blunting of the day-night 
variations of the RR interval 1911 and LF and HF spectral 
components 191,921 present in a period ranging from days 
to a few weeks after the acute event. In post MI patients 
with a very depressed HRV, most of the residual energy 
is distributed in the VLF frequency range below 0 03 Hz, 
with only a small respiration-related HF 1931 . These char- 
acteristics of the ^spectral profile are similar to those 
observed in advanced cardiac failure or after cardiac 
transplant, and are likely to reflect either diminished 
responsiveness of the target organ to neural modulatory 
inputs 1 * 21 or a saturating influence on the sinus node of a 
persistently high sympathetic tone 1281 . 

Diabetic neuropathy 

In neuropathy associated with diabetes mellitus, charac- 
terized by alteration of small nerve fibres, a reduction in 
time-domain parameters of HRV seems not only to 
carry negative prognostic value but also to precede the 
clinical expression of autonomic neuropathy 194 " 971 . In 
diabetic patients without evidence of autonomic neuro- 
pathy, reduction of the absolute power of LF and HF 
during controlled conditions was also reported* 961 . How- 
ever, when the LF/HF ratio was considered or when LF 
and HF were analysed in normalized units, no signifi- 
cant difference in comparison to normals was present 
Thus, the initial manifestation of this neuropathy is 
likely to involve both efferent limbs of the autonomic 
nervous system 196,981 . 

Cardiac transplantation 

A very reduced HRV with no definite spectral compo- 
nents was reported in patients with a recent heart 
transplant 197, * 100i . The appearance of discrete spectral 
components in a few patients is considered to reflect 
cardiac re-innervauW 1011 . This re-innervation may 
occur as early as 1 to 2 years post transplantation 
and is usually of sympathetic origin. Indeed, the corre- 
lation between the respiratory rate and the HF compo- 
nent of HRV observed in some transplanted patients 
indicates that a non-neural mechanism may also contrib- 
ute to generate respiration-related rhythmic oscilla- 
tion 11001 . The initial observation of identifying patients 
developing an allograft rejection according to changes in 
HRV could be of clinical interest but needs further 
confirmation. 

Myocardial dysfunction 

A reduced HRV has been consistently observed in 
patients with cardiac failure* 24 .WUW-1MT in this condi- 
tion characterized by signs of sympathetic activation. 



such as faster heart rates and high levels of circulating 
cathecolamines, a relationship between changes in HRV 
and the extent of left ventricular dysfunction was con- 
troversially reported 1102,1 MJ . In fact, whereas the reduc- 
tion in time domain measures of HRV seemed to parallel 
the severity of the disease, the relationship between 
spectral components and indices of ventricular dysfunc- 
tion appears to be more complex. In particular, In most 
patients with a very advanced phase of the disease and 
with a drastic reduction in HRV, a LF component could 
not be detected despite clinical signs of sympathetic 
activation. Thus, in conditions characterized by marked 
and unopposed persistent sympathetic excitation, the 
sinus node seems to drastically diminish its responsive- 
ness to neural inputs 11041 . 

Tetraplegia 

Patients with chronic complete high cervical spinal cord 
lesions have intact efferent vagal and sympathetic neural 
pathways directed to the sinus node. However, spinal 
sympathetic neurons are deprived of modulatory control 
and in particular of baroreflex supraspinal inhibitory 
inputs. For this reason, these patients represent a unique 
clinical model with which to evaluate the contribution of 
supraspinal mechanisms in determining the sympathetic 
activity responsible for low frequency oscillations of 
HRV. It has been reported 11071 that no LF could be 
detected in tetraplegic patients, thus suggesting the criti- 
cal role of supraspinal mechanisms in detennining the 
01 Hz rhythm. Two recent studies, however, have indi- 
cated that an LF component can also be detected in 
HRV and arterial pressure variabilities of some tetra- 
plegic patients 11081091 . While Koh et alS m attributed the 
LF component of HRV to vagal modulations, Guzzetti 
et a/. 11091 attributed the same component to sympathetic 
activity because of the delay with which the LF com- 
ponent appeared after spinal section, suggesting an 
emerging spinal rhythmicity capable of modulating 
sympathetic discharge. 

Modifications of HRV by specific 
interventions 

The rationale for trying to modify HRV after MI stems 
from multiple observations indicating that cardiac 
mortality is higher among post-MI patients who have a 
more depressed HRV 193110 *. The inference is that inter- 
ventions that augment HRV may be protective against 
cardiac mortality and sudden cardiac death. Although 
the rationale for changing HRV is sound, it contains 
also the inherent danger of leading to the unwarranted 
assumption that modification of HRV translates directly 
into cardiac protection, which may not be the case 111 ". 
The target is the improvement of cardiac .electrical 
stability, and HRV is just a marker of autonomic 
activity. Despite the growing consensus that increases in 
vagal activity can be beneficial 11 1Z1 , it is not yet known 
how much vagal activity (or its markers) has to increase 
m order to provide adequate protection. 
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Beta-adrenergic blockade and HRV 
The data on the effect of ^-blockers on HRV in post- 
Mi patients are surprisingly scanty 11131141 . Despite 
the observation of statistically significant increases, the 
actual changes are very modest. However, it is of note 
that /^-blockade prevents the rise in the LF component 
observed in the morning hours 11141 . In conscious post-MI 
dogs, ^-blockers do not modify HRV 1 " 51 . The un- 
expected observation that, prior to MI, /^-blockade 
increases HRV only in the animals destined to be at low 
risk for lethal arrhythmias post-MI 11151 may suggest 
novel approaches to post-MI risk stratification. 

Antiarrhythmic drugs and HRV 
Data exist for several antiarrhythmic drugs. Flecainide 
and propafenone, but not amiodarone, were reported to 
decrease time-domain measures of HRV in patients with 
chronic ventricular arrhythmia 11161 . In another study 11171 
propafenone reduced HRV and decreased LF much 
more than HF, resulting in a significantly smaller 
LF/HF ratio. A larger study" 181 confirmed that flecain- 
ide, and also encamide and moricizine, decreased HRV 
in post-MI patients but found no correlation between 
the change in HRV and mortality during follow-up. 
Thus, some antiarrhythmic drugs associated with 
increased mortality can reduce HRV. However, it is not 
known whether these changes in HRV have any direct 
prognostic significance. 

Scopolamine and HRV 

Low dose muscarinic receptor blockers, such as atropine 
and scopolamine, may produce a paradoxical increase in 
vagal efferent activity, as suggested by a decrease in 
heart rate. Different studies examined the effects of 
transdermal scopolamine on indices of vagal activity in 
patients with a recent MI ,11& ~ 1221 and with congestive 
heart failure 11231 . Scopolamine markedly increases HRV, 
which indicates that pharmacological modulation of 
neural activity with scopolamine may effectively increase 
vagal activity. However, efficacy during long-term treat- 
ment has not been assessed. Furthermore, low dose 
scopolamine does not prevent ventricular fibrillation due 
to acute myocardial ischaemia in post-MI dogs 11241 . 

Thrombolysis and HRV 

The effect of thrombolysis on HRV (assessed by 
pNN50), was reported in 95 patients with acute MI 1125 * 
HRV was higher 90 rnin after thrombolysis in the 
patients with patency of the infarct-related artery. How- 
ever, this difference was no longer evident when the 
entire 24 h were analysed. 

Exercise training and HRV 

Exercise training may decrease cardiovascular mortality 
and sudden cardiac death 11261 . Regular exercise training 
is also thought capable of modifying the autonomic 
balance 1127,12 *. A recent experimental study, designed to 
assess the effects of exercise training on markers of vagal 
activity, has simultaneously provided information on 
changes in cardiac electrical stability 11291 . Conscious 



dogs, documented to be at high risk by the previous 
occurrence of ventricular fibrillation during acute myo- 
cardial ischaemia, were randomly assigned to 6 weeks of 
either daily exercise training or cage rest followed by 
exercise training 11291 . After training, HRV (SDNN) in- 
creased by 74% and all animals survived a new ischaemic 
test. Exercise training can also accelerate recovery of the 
physiological sympatho-vagal interaction, as shown in 
post-MI patients"* 01 . 

Clinical use of heart rate variability 

Although HRV has been the subject of numerous clini- 
cal studies investigating a wide spectrum of cardiological 
and non-cardiological diseases and clinical conditions, a 
general consensus of the practical use of HRV in adult 
medicine has been reached only in two clinical scenarios. 
Depressed HRV can be used as a predictor of risk after 
acute MI and as an early warning sign of diabetic 
neuropathy. 

Assessment of risk after acute myocardial 
infarction 

The observation 1121 that in patients with an acute MI the 
absence of respiratory sinus arrhythmias is associated 
with an increase in 'in-hospital' mortality represents the 
first of a large number of reports 116,931311 which have 
demonstrated the prognostic value of assessing HRV to 
identify high risk patients. 

Depressed HRV is a powerful predictor of mor- 
tality and of arrhythmic complications (e.g. sympto- 
matic sustained ventricular tachycardia) in patients 
following acute MI 116,1311 (Fig. 8). The predictive value 
of HRV is independent of other factors established for 
post-infarction risk stratification, such as depressed 
left ventricular ejection fraction, increased ventricular 
ectopic activity, and presence of late potentials. For 
prediction of all-cause mortality, the value of HRV is 
similar to that of left ventricular ejection fraction, but 
HRV is superior to left ventricular ejection fraction in 
predicting arrhythmic events (sudden cardiac death and 
ventricular tachycardia) 11313 . This permits speculation 
that HRV is a stronger predictor of arrhythmic mortal- 
ity rather than non-arrhythmic mortality. However, 
clear differences between HRV in patients suffering from 
sudden and non-sudden cardiac death after acute MI 
have not been observed. Nevertheless, this might also be 
related to the nature of the presently used definition of 
sudden cardiac death 11321 , which is bound to include not 
only patients suffering from arrhythmia-related death 
but also fatal reinfarctions and other cardiovascular 
events. 

The value of both conventional time-domain 
and frequency-domain parameters have been fully 
assessed in several independent prospective studies, but 
because of using optimized cut-off values defining nor- 
mal and depressed HRV, these studies may slightly 
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Figure 6* Cumulative survival of patients after ML (a) Shows 
survival of patients stratified according to 24-h SDNN values in 
three groups with cut-off points of 50 and 100 ms. (Reprinted with 
permission* 161 , (b) Shows similar survival curves of patients strati- 
fied according to 24-h HRV triangular index values with cut-off 
points of 15 and 20 units, respectively. (Data of St. George's 
Post-infarction Research Survey Programme.) 



overestimate the predictive value of HRV. Nevertheless, 
the confidence intervals of such cut-off values are rather 
narrow because of the sizes of the investigated popula- 
tions. Thus, the observed cut-off values of 24-h measures 
of HRV, e.g. SDNN <50 ms and HRV triangular index 
<15 for highly depressed HRV, or SDNN < 100 ms and 
HRV triangular index <20 for moderately depressed 
HRV are likely to be broadly applicable. 

It is not known whether different indices of HRV 
(e.g. assessments of the short- and long-term compo- 
nents) can be combined in a multivariate fashion in 
order to improve post-infarction risk stratification. 
There is a general consensus, however, that combination 
of other measures with the assessment of overall 24-h 
HRV is probably redundant 

Pathophysiological considerations 

It has not been established whether depressed HRV is 

part of the mechanism of increased post-infarction mor- 



tality or is merely a marker of poor prognosis. The data 
suggest that depressed HRV is not a simple reflection of 
the sympathetic overdrive and/or vagal withdrawal due 
to poor ventricular performance but that it also reflects 
depressed vagal activity which has a strong association 
with the pathogenesis of ventricular arrhythmias and 
sudden cardiac death 1113 . 

Assessment of HRV for risk stratification after acute 
myocardial infarction 

Traditionally, HRV used for risk stratification after MI 
has been assessed from 24-h recordings. HRV measured 
from short-term electrocardiogram recordings also 
provides prognostic information for risk stratification 
following MI but whether it is as powerful as that from 
24-h recordings is uncertam* 133 -™ 51 . HRV measured 
from short-term recordings is depressed in patients at 
higji risk; the predictive value of depressed HRV in- 
creases with increasing the length of recording. Thus, the 
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Figure 9 Comparison of positive predictive characteris- 
tics ofHRV (solid lines) and of combinations of HRV with 
left ventricular ejection fraction (dashed lines) and of 
HRV with left ventricular ejection fraction and ectopic 
counts on 24-h ECGs (dotted lines) used for identification 
of patients at risk of 1-year cardiac mortality (a) and 
1-year arrhythmic events (sudden death and/or sympto- 
matic sustained ventricular tachycardia (b) after acute 
myocardial infarction. (Data of St George's Post- 
infarction Research Survey Programme.) 



use of nominal 24-h recordings may be recommended 
for risk stratification studies after MI. On the other 
hand, the assessment of HRV from short-term record- 
ings can be used for initial screening of survivors of 
acute MI Il3eI . Such an assessment has similar sensitivity 
but lower specificity for predicting patients at high risk 
compared to 24-h HRV. 

Spectral analysis of HRV in survivors of MI 
suggested that the ULF and VLF components carry the 
highest predictive value* 93 *. As the physiological correlate 
of these components is unknown and as these compo- 
nents correspond to up to 95% of the total power which 
can be easily assessed in the time-domain, the use of 
individual spectral components of HRV for risk stratifi- 
cation after MI is not more powerful than the use of 
those time-domain methods which assess overall HRV. 

Development of HRV after acute myocardial infarction 
The time after acute MI at which the depressed HRV 
reaches the highest predictive value has not been inves- 
tigated comprehensively. Nevertheless, the general con- 
sensus is that HRV should be assessed shortly prior to 
hospital discharge, i.e. approximately 1 week after index 
infarction. Such a recommendation also fits well into the 
common practice of hospital management of survivors 
of acute ML 



Heart rate variability is decreased early after 
acute MI and begins to recover within a few weeks; it is 
maximally but not fully recovered by 6 to 12 months 
after MP 11371 . Assessment of heart rate variability at 
both the early stage of MI (2 to 3 days after acute MI)* 841 
and pre-discharge from hospital (1 to 3 weeks after acute 
MI) offers important prognostic information. Heart rate 
variability measured late (1 year) after acute MI also 
predicts further mortality 11 381 . Data from animal models 
suggest that the speed of HRV recovery after MI 
correlates with subsequent risk 11151 . 

HRV used for multivariate risk stratification 
The predictive value of heart rate variability alone Is 
modest, but combination with other techniques substan- 
tially improves its positive predictive accuracy over a 
clinically important range of sensitivity (25% to 75%) for 
cardiac mortality and arrhythmic events (Fig. 9). 

Improvements in the positive predictive accuracy 
over the range of sensitivities have been reported 
for combinations of HRV with mean heart rate, left 
ventricular ejection fraction, frequency of ventricular 
ectopic activity, parameters of high resolution electro- 
cardiograms (e.g. presence or absence of late potentials), 
and clinical assessment 11331 . However, it is not known 
which other stratification factors are the most practical 
and most feasible to be combined with HRV for 
multifactorial risk stratification. 

Systematic multivariate studies of post MI risk 
stratification are needed before a consensus can be 
reached and before a combination of HRV with other 
variables of proven prognostic importance can be rec- 
ommended. Many aspects that are not relevant for 
univariate risk stratification need to be examined: it is 
not obvious whether the optimum cut-off values of 
individual risk factors known from univariate studies are 
appropriate in a multivariate setting. Different multi- 
variate combinations are probably needed for optimiz- 
ing predictive accuracy at different ranges of sensitivity. 
Stepwise strategies should be examined to identify 
optimum sequences of performing individual tests used 
in multivariate stratification. 

Summary and recommendations for interpreting 
predictive value of depressed HRV after acute 
myocardial infarction 

The following facts should be noted when exploiting 
HRV assessment in clinical studies and/or trials 
involving survivors of acute myocardial infarction. 

Depressed HRV is a predictor of mortality and 
arrhythmic complications independent of other 
recognised risk factors. 

There is a general consensus that HRV should be 
measured approximately 1 week after index infarction. 

Although HRV assessed from short-term record- 
ings provides prognostic Information, HRV measured in 
nominal 24-h recordings is a stronger risk predictor. 
HRV assessed from short-term recordings may be used 
for initial screening of all survivors of an acute MI. 
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No currently recognised HRV measure provides 
better prognostic information than the time-domain 
HRV measures assessing overall HRV (e.g. SDNN or 
HRV triangular index). Other measures, e.g. ULF of 
entire 24-h spectral analysis, perform equally well. A 
high risk group may be selected by the dichotomy limits 
of SDNN <50 ms or HRV triangular index < 15. 

For clinically meaningful ranges of sensitivity, 
the predictive value of HRV alone is modest, although it 
is higher than that of any other so far recognised risk 
factor. To improve the predictive value, HRV may be 
combined with other factors. However, optimum set of 
risk factors and corresponding dichotomy limits have 
not yet been established.. 



Assessm*At of diabetic neuropathy 

As a complication of diabetes mellitus, autonomic 
neuropathy is characterized by earry and widespread 
neuronal degeneration of small nerve fibres of both 
sympathetic and parasympathetic tracts 11401 . Its clinical 
manifestations are ubiquitous with functional impair- 
ment and include postural hypotension, persistent tachy- 
cardia, gustatory sweating, gastroparesis. bladder atony 
and nocturnal diarrhoea. Once clinical manifestations of 
diabetic autonomic neuropathy (DAN) supervene, the 
estimated 5-year mortality is approximately 50% 11411 . 
Thus, early subclinical detection of autonomic dysfunc- 
tion is important for risk stratification and subsequent 
management Analyses of short-term and/or long-term 
HRV have proven useful in detecting DAN 1961 42-1 471 . 

For the patient presenting with a real or suspect 
DAN there are three HRV methods from which to 
choose: (a) simple bedside RR interval methods, (b) 
long-term time-domain measures, which are more 
sensitive and more reproducible than the short-term 
tests, and (c) frequency-domain analysis performed 
under short-term steady state conditions and which are 
useful in separating sympathetic from parasympathetic 
abnormalities. 

Long-term time domain measures 
HRV computed from 24-h Holter records are more 
sensitive than simple bedside tests (e.g. Valsava 
manoeuvre, orthostatic test, and deep breathing 1111 ) for 
detecting DAN. Most experience has been obtained with 
the NN50 D441 and SDSD (see Table l) 0451 methods. 
Using the NN50 count, where the lower 95% confidence 
interval for total counts range from 500 to 2000 depend- 
ing on the age. about half of diabetic patients will 
demonstrate abnormally low counts per 24 h. Moreover, 
there is a strong correlation between the percentage 
of patients with abnormal counts and the extent of 
autonomic neuropathy determined from conventional 
measures, 

Besides their increased sensitivity, these 24-h 
time domain methods are strongly correlated with other 
established HRV measurements and have been found 
to be reproducible and stable over time. Similar to 



survivors of MI, patients with DAN are also predis- 
posed to poor outcomes such as sudden death but it 
remains to be determined whether the HRV measures 
confer prognostic information among diabetics. 

Frequency domain measures 

The following abnormalities in frequency HRV analy- 
sis are associated with DAN: (a) reduced power in 
all spectral bands which is the most common find- 
^{95.146-1481 (5) failure to increase LF on standing, 
which is a reflection of impaired sympathetic response or 
depressed baroreceptor sensitivity 196,1471 ; (c) abnormally 
reduced total power with unchanged LF/HF ratio 196 *, 
and (d) a leftward shift in the LF central frequency, 
the physiological meaning of which needs further 
elucidation 1147 ". 

In advanced neuropathic states, the resting 
supine power spectrum often reveals extremely low 
amplitudes of all spectral components making it difficult 
to separate signal from noise 196 * 146 1471 . It is therefore 
recommended that an intervention such as standing or 
tilt be Included. Another method to overcome the low 
signal to noise ratio is to introduce a coherence function 
which utilizes the total power coherent with one or the 
other frequency band 11461 . 



Other clinical potential 

Selected studies investigating HRV in other cardio- 
logical diseases are listed in Table 4. 



Future possibilities 

Development of HRV measurement 

The currently available time-domain methods predomi- 
nantly used to assess the long-term profile of HRV are 
probably sufficient for this purpose. Improvements are 
possible, especially in terms of numerical robustness. 
The contemporary non-parametric and parametric 
spectral methods are probably sufficient to analyse 
short-term ECGs without transient changes of heart 
period modulations. 

Apart from the need to develop numerically 
robust techniques suitable for fully automatic measure- 
ment (the geometrical methods are only one possibility 
in this direction), the foDowing three areas deserve 
attention. 

Dynamics and transient changes of HRV 
The present possibilities of characterizing the quantify- 
ing the dynamics of the RR interval sequence and 
transient changes of HRV are sparse and still under 
mathematical development. However, it ' may be 
assumed that proper assessment of HRV dynamics will 
lead to substantial improvements in our understanding 
of both the modulations of heart period and their 
physiological and pathophysiological correlates. 
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Table 4 Summary of selected 
myocardial infarction 



studies investigating the clinical value of HRV In cardiological diseases other than 



Disease state 



Author 
of study 



Population 
(no. of patients} 



Investigation 
parameter 



Clinical finding 



Potential value 



Hypertension Guzzetti 
1991'" 31 



Congestive heart * Saul 

1988i |5SJ 



CasoJo 
1989 liozj 



Binkley 
1991 tift 



Kienzte 
1992^ 



Townend 
1992l ,S3 l 

Binkley 
I993' 1 *" 



Woo 
1994 n551 



49 hypertensive Spectral AR 
30 normals 



41 borderline 
hypertensive 
34 hypertensive 
54 normals 

25 chronic CHF 
NYHA III. IV 
21 normals 



Spectral FFT 



t LF found in 
hypertensives as compared 
to normals with blunting 
of drcadian patters 

Reduced parasympathetic 
in hypertensive patients 



Spectral i spectral power all 

Blackman-Turkey frequencies, especially 

15min >0-04HzinCHF 
acquisition 



20 CHF Time-domain RR Low HRV 

NYHA II. Ill, IV interval histogram 
20 normals with 24 h-Holter 



10 dilated 
cardiomyopathy 
(EF 14 to 40%) 
10 normals 



Spectral FFT 
4 min supine 
acqui si ti o n 



23 CHF Spectral FFT 

NYHA II. III. IV Time-domain 
24-48 h-Holter 



12 CHF 
NYHA HI. IV 



Time domain 
24 h-Holter 



13 CHF Spectral FFT 

NYHA II. Ill 4 min supine 
acquisition 



21 CHF 
NYHA m 



Polncare plots 
Time-domain 
24 h-Holter 



I m high frequency power 
(>0-l Hz) in CHF 
T LF/HF 



Alterations of HRV not 
tightly linked to severity of 
CHF 

i HRV was related to 
sympathetic excitation 

HRV T during ACE 
inhibitor treatment 

12 weeks of ACE inhibitor 
treatment t high 
frequency HRV 

Complex plots are 
associated with t 
norepinephrine levels and 
greater sympathetic 
activation 



Hypertension is 
characterized by a 
depressed drcadian 
rhythmlcity of LF 

Support the use of 
non-pathological therapy of 
hypertension mat | vagal 
tone (e.g. exercise) 

In CHF. there is i vagal, but 
relatively preserved 
sympathetic modulation of 
HR 

Reduced vagal activity in 
CHF patients 

Withdrawal of 
parasympathetic tone 
observed in CHF 
CHF has imbalance of 
autonomic tone with X 
parasympathetic and a 
predominance of 
sympathetic tone 



A significant augmentation 
of parasympathetic tone 
was associated with ACE 
inhibitor therapy 

Polncare plots may assist 
analysis of sympathetic 
influences 



Heart Alexopoulos 19 transplant 

transplantation 198$** 10 i 



Sands 



17 transplant 
6 normals 



Time-domain 
24 h-Holter 



Spectral FFT 
15 min supine 
acquisition 



Chronic mitral Stein 38 chronic mitral Spectral FFT 

regurgitation 1993 11571 regurgitation Time-domain 

24 h-Holter 



Mitral valve 
prolapse 



Marangoni 
1993 |,5 « 



39 female mitral 
valve prolapse 
24 female 
controls 



Spectral AR 
10 min supine 
acquisition 



Reduced HRV in 
denervated donor hearts; 
recipient innervated 
hearts had more HRV 

HRV from 002 to 1-0 Hz- 
90% reduced 



HR and measures of 
ultrakw frequency by 
SDANN correlated with 
ventricular performance 
and predicted clinical events 

MVP patients had 1 high 
frequency 



Patients with rejection 
documented biopsy show 
significantly more variability 

May be prognostic indicator 
of atrial fibrillation, mortality, 
and progression to valve 
surgery 



MVP patients had low 
tone 



Continued 
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Table 4 Continued. 



Disease state 



Author Population 
of study {no. of patients) 



Investigation 
parameter 



Clinical finding 



Potential value 



Cardiomyopathies Courrihan 



Sudden death or 
cardiac arrest 



Dougherty 16 CA survivors Spectral AR 
1992^ 60! 5 CA 

nonsurvivors 

5 normals 

HuikurT 22 CA survivors 
1992 ,,€11 22 control 



Algra 
I9931 1 "* 



Ventricular 
arrhythmias 



Myers 
1986 l,6Z » 



Martin 
1988 13631 



Vybiral 



104 HCM Spectral FFT Global and specific vagal 

Time-domain tone measurements of 
24 h-Holter HRV were j. in 

symptomatic patients 

HRV as measured by low 
Time-domain frequency power and 
24 h-Holter SDNN were significantly 

related to 1 year mortality 

Spectral AR J High frequency power in 

Time-domain CA survivors — low 
24 h-Holter frequency power did not 

discriminate CA survivors 
Circadian pattern of HRV 
found m all patients 

Tune-domain I short-term variation 
24 h-Holter (a05-O-50 Hz) 

independently increased 
the risk of SD by a factor 
2-6 

4 long-terra variation 
(0 02-0-05 Hz) increased 
the risk of SD by a factor 
of2 

Both time and frequency 
domain indices separated 
normals from SD patients 
I HF power (0-35-0-5 Hz) 
was the best separator 
between heart disease 
patients with and without SD 

SDNN index significantly 
tower in SD patients 



HRV does not add to the 
predictive accuracy of 
known risk factors in HCM 



HRV is clinically useful to 
risk stratify CA survivors for 
1 year mortality 



193 SD cases 
230 

symptomatic 
patients 



6 normals 
12 patients with 
structural heart 
disease (6 with 
and 6 without 
SD) 



HRV may be used to 
estimate the risk of SD 



Time and 
frequency domain 
24 h-Holter 



20 normals 
5 patients 
experiencing SD 
during Hotter 
monitoring 

24 VF 
19 CAD 



Time-domain 
24 h-Holter 



Hulkuri 18 VT or CA 
1993« 1BS > 



Hohnloser 14 post MI with Spectral FFT 
1994 [1561 VF or sustained Time-domain 
14 post MX 24 h-Holter 



Supraventricular Kocovic 
arrhythmias 1993 ,1CT1 



HF power may be a useful 
predictor of SD 



Time domain indices may 
Identify increased risk of SD 



Time-domain HRV Indices do not 
24 h-Holter change consistently 

before VF 

Spectral AR all power spectra of HRV 

24 h-Holter were significantly J. before 

the onset of sustained VT 
than before nonsustained VT 

HRV of post MI CA 
survivors do not differ 
from other post MI 
patients they differ strikingly 
in term of bar o reflex 
sensitivity 



64SVT Spectral FFT t HR. i HRV and I 

Time-domain parasympathetic 

acquisition components after 

24 h-Holter radioftequency ablation 



A temporal relation exists 
between the decrease of 
HRV and the onset of 
sustained VT 

Baroreflex sensitivity, not 
HRV. distinguished post MI 
patients with and without VF 
and VT 



Parasympathetic ganglia 
and fibres may be more 
dense in the mid and 
anterior low septum 



AR « autnregressi ve; C A -cardiac arrest: C AD =■ coronary artery disease; CHF= congestive heart failure; EF-ejection fraction; FFT = fast 
Fourier transform; HCM = hypertrophic cardiomyopathy; HF=high frequency; HRV = heart rate variability; LF=low frequency; 
NYHA=New York Heart Association classification; SD -sudden death: SVT=supraventricular tachycardia; VF=ventricular fibrillation; 
VT- ventricular tachycardia. 
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It remains to be seen whether the methods of 
non-linear dynamics will be appropriate for the 
measurement of transient changes in RR intervals or 
whether new mathematical models and algorithmic 
concepts will be needed to tailor the principles of 
measurement more closely to the physiological nature of 
cardiac period ograms. In any case, the task of assessing 
transient changes in HRV seems to be more important 
than further refinements of the current technology used 
to analyse stable stages of heart period modulations. 

PP and RR intervals 

Little is known about the interplay between the PP 
and PR autonomic modulations. For these reasons, the 
sequence of PP~ intervals should also be studied 11681 . 
Unfortunately, precise location of a P wave fiducial 
point is. almost impossible to achieve in surface ECGs 
recorded \ with the current technology. However, 
developments in the technology should allow PP interval 
and PR interval variability to be investigate in future 
studies. 

Multisignal analysis 

The modulations of heart periods are naturally not 
the only manifestation of the autonomic regulatory 
mechanisms. Currently, commercial or semi-commercial 
equipment exists which enables simultaneous recording 
of ECG, respiration, blood pressure, etc. However, in 
spite of the ease with which the signals can be recorded, 
no widely accepted method exists for comprehensive 
multisignal analysis. Each signal can be analysed separ- 
ately, e.g. with parametric spectral methods, and the 
results of the analysis compared. Analysis of coupling 
between physiological signals allow the properties of the 
coupling to be measured* 1 69 ~ 1741 . 

Studies needed to improve physiological 
understanding 

Efforts should be made to find the physiological corre- 
lates and the biological relevance of various HRV 
measures currently employed. In some cases, e,g. the HF 
component, this has been achieved. In other cases, 
e.g. the VLF and ULF components, the physiological 
correlates are still largely unknown. 

This uncertain knowledge limits the interpreta- 
tion of associations between these variables and the risk 
of cardiac events. The use of markers of autonomic 
activity is very attractive. However, unless a tenable 
mechanistic link between these variables and cardiac 
events is found, there is an inherent danger of concen- 
trating therapeutic efforts on the modification of these 
markers Iuuf2J . This may lead to incorrect assumptions 
and serious misinterpretations. 

Possibilities of future clinical utility 
Normal standards 

Large prospective population studies with longitudinal 
follow-up are needed to establish normal HRV stan- 



dards for various age and gender subsets* 1101 . Recently, 
investigators from the Framingham Heart Study 
reported on the time- and freq uency-dornain measures 
of HRV in 736 elderly subjects, and the relationship of 
these HRV measures to all-cause mortality during 4 
years of foHow-up [1751 . These investigators concluded 
that HRV offers prognostic information independent of 
and beyond that provided by traditional risk factors. 
Additional population-based HRV studies involving the 
full age spectrum in males and females need to be 
performed. 

Physiologic phenomena 

It would be of interest to evaluate HRV in various 
arcadian patterns such as normal day-night cycles, 
sustained reversed day-night cycles (evening-night shift 
work), and transiently altered day-night cycles, such as 
might occur with international travel. The autonomic 
fluctuations occurring during various stages of sleep 
including rapid eye movement (REM) sleep have been 
studied in only a few subjects. In normal subjects, the 
HF vagal component of the power spectrum is aug- 
mented only during non-REM sleep, whereas in post-MI 
patients with this increase in HF is absent 1176 *. 

The autonomic nervous system response to ath- 
letic training and rehabilitative exercise programmes 
after various disease states is thought to be a condition- 
ing phenomenon. HRV data should be useful In under- 
standing the chronological aspects of training and the 
time to optimal conditioning as it relates to the auto- 
nomic influences on the heart. Also, HRV may provide 
important information about deconditioning with pro- 
longed bed rest, with weightlessness and with the zero g 
that accompanies space flight. 

Pharmacological responses 

Many medications act directly or indirectly on the 
autonomic nervous system, and HRV can be used to 
explore the influence of various agents on sympathetic 
and parasympathetic activity. It is known that para- 
sympathetic blockade with full dose atropine produces 
marked diminution of HRV. Low dose scopolamine has 
vagotonic influences and is associated with Increased 
HRV, especially in the HR range, ^-adrenergic blockade 
was observed into increase HRV and to reduce the 
normalized units of the LF component? 151 . Considerably 
more research is needed to understand the effects and 
clinical relevance of altered vagotonic and adrenergic 
tone on total HRV power and its various components in 
health and disease. 

At present, few data exist on the effects of 
calcium channel blockers, sedatives, anxiolytics, analge- 
sics, anaesthetics, antiarrhythmic agents, narcotics, 
and chemotherapeutic agents such as vincristine and 
doxorubicin on HRV. 

Risk stratification 

Both time and frequency measures of HRV calculated 
from long 24-h and short 2 to 15-min ECG recordings 
have been used to predict time to death after MI, as well 
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as the risk of all-cause mortality and sudden caitijac 
death In patients with structural heart disease 1 ■ ' 
and a number of other pathophysiological condi- 
tions* 1771 . Using diagnostic instruments that can measure 
HRV. together with the frequency and complexity 
of ventricular arrhythmias, signal-averaged ECG. ST 
segment variability, and repolarization heterogeneity, it 
should be possible to markedly improve the identifica- 
tion of patients at risk for sudden cardiac death and 
arrhythmic events. Prospective studies are needed to 
evaluate the sensitivity, specificity, and predictive 
accuracy of combined testing. 

Fetal and neonatal heart rate variability is an 
important area of investigation, and it might provide 
early information about fetal and neonatal distress 
and identify those at risk for the sudden infant death 
syndrome. Mostiofthe preliminary work in this field was 
carried out in the -early 1980s before the more sophisti- 
cated power spectral techniques became available. In- 
sight into autonomic maturation in the developing fetus 
might also be possible through the proper application of 
these techniques. 

Disease mechanisms 

A fertile area of research is to use HRV techniques to 
explore the role of autonomic nervous system alterations 
in disease mechanisms, especially those conditions in 
which sympathovagal factors are thought to play an 
important role. Recent work suggests that alterations in 
autonomic innervation to the developing heart might 
be responsible for some forms of the long QT syn- 
drome" 781 . Fetal HRV studies in pregnant mothers with 
this disorder is certainly feasible and might be very 
informativef 1791 . 

The role of the autonomic nervous system in 
essential hypertension is an important area of investiga- 
tion 11801 . The question regarding the primary or second- 
ary role of enhanced sympathetic activity in essential 
hypertension might be answered by longitudinal studies 
of subjects who are initially normotensive. Does essen- 
tial hypertension result from augmented sympathetic 
activity with altered responsiveness of neural regulatory 
mechanisms? 

Several primary neurological disorders including 
Parkinson's disease, multiple sclerosis. Guillain-Barre 
syndrome, and orthostatic hypotension of the Shy- 
Drager type are associated with altered autonomic func- 
tion. In some of these disorders, changes in HRV may be 
an early manifestation of the condition and may be 
useful in quantltating the rate of disease progression 
and/or the efficacy of therapeutic interventions. This 
same approach may also be useful in the evaluation of 
secondary autonomic neurologic disorders that accom- 
pany diabetes mellitus, alcoholism, and spinal cord 
injuries. 

Conclusion 

Heart rate variability has considerable potential to as- 
sess the role of autonomic nervous system fluctuations in 
normal healthy individuals and in patients with various 



cardiovascular and non-cardiovascular disorders. HRV 
studies should enhance our understanding of physiologi- 
cal phenomena, the actions of medications, and disease 
mechanisms. Large prospective longitudinal studies are 
needed to determine the sensitivity, specificity, and pre- 
;tive value of HRV in the Identification of individuals 
risk for subsequent morbid and mortal events. 
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Appendix A 

Normal values of standard measures of heart 
rate variability 

As no comprehensive investigations of all HRV Indices 
in large normal populations have yet been performed. 
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Variable 



Units 



Normal values 
(mean±SD) 
Time domain analysis of nominal 24 h 081 ' 



SDNN 


ms 


141 ±39 


SDANN 


ms 


127 ±35 


RMSSD 


ms 


27 ±12 


HRV triangular index 




37 ±15 


Spectral analysis of stationary supine 5-min recording 


total power 


ms 2 


3466 ±1018 


LF 


ms 2 


1170±416 


HF 


ms 2 


975 ±203 


LF 


n.u. 


54±4 


HF 


mi. 


29±3 


LF/HF ratio 

: » 




1.5-20 



some of the normal values listed in the following table 
were obtained from studies involving small number of 
subjects. The values should therefore have been consid- 
ered as approximate and no definite clinical conclusions 
should be based on them. The adjustment of normal 
limits for age, sex, and environment which is also needed 
has been omitted here because of the limited sources of 



The table lists only values of those measures of 
HRV which might be suggested for standardisation of 
further physiological and clinical studies. 



Appendix B 

Suggestion of procedures for testing of 
commercial equipment designed to measure 
heart rate variability 

Concept 

In order to achieve comparable accuracy of measure- 
ments reported by different commercial equipment, each 
device should be tested independently of the manufac- 
turer (e.g; by an academic institution). Each test should 
involve several short-term and, if applicable, long-term 
test recordings with precisely known HRV parameters 
and with different morphological characteristics of the 
ECG signaL If the involvement of the manufacturer is 
required during the testing procedure (e.g. for manual 
editing of the labels of QRS complexes), the manufac- 
turer must be blinded in respect of both the true HRV 
parameters of the testing recordings and the features 
used to obtain the signal. In particular, when the results 
of the test are disclosed to the manufacturer for further 
improvement of the device or other wise, new tests 
should involve a completely new set of test recordings. 

Techoicai requirements 

Each device should be tested comprehensively including 
all its parts. In particular, the test should involve both 



the recording and the analytical part of the device. An 
appropriate technology should be used to record a fully 
reproducible signal with precisely known HRV param- 
eters, e.g. the test signal should be computer and/or 
hardware generated. Both brand new recorders as well 
as recorders which have been routinely used and rou- 
tinely serviced for approximately a half of their lifetime 
should be used in the tests if this is feasible (the testing 
should not be delayed for newly introduced systems). 
If a manufacturer claims that the device is capable 
of analysing ECG records (such as Holter tapes) ob- 
tained with recorders of other manufacturers, each 
combination should be tested Independently. 

As the analysis of HRV by implantable devices 
may be foreseen, similar procedures as described further 
should be used to generate simulated intracardiac sig- 
nals. If feasible, implantable devices with fully charged 
batteries as well as devices with partly discharged 
batteries should be tested. 

Test recordings 

It is intrinsically difficult to know precisely the HRV 
parameters of any real ECG recordings independently of 
equipment used to analyse the recording. Therefore, 
simulated ECG signals are preferable. However, the 
morphology of such simulated ECG signals as well as 
the HRV characteristics must closely reflect the mor- 
phology of real recordings. The discrete frequency used 
to generate such signals must be substantially higher 
than the sampling frequency of the tested device. Fea- 
tures which should be introduced into such recordings 
should include different factors known to influence or 
potentially influence the precision of HRV assessment, 
e.g. variable noise levels, variable morphology of QRS 
signals which may cause jitter of the fiducial point, 
randomly alternating noise in different channels of the 
signal, gradual and abrupt changes of HRV character- 
istics, and different frequencies of atrial and ventricular 
ectopic beats with realistic morphologies of the signal. 

The quality of records on magnetic tape based 
systems may not be constant during long-term recording 
due to spool torque control, back tension, and other 
factors. Performance of all recorders can be influenced 
by changes of the outside environment. Long-term 
(e.g. full 24-h test) rather than short-term tests should 
therefore be used. 

Testing procedures 

Each device or each configuration of the device should 
be tested with several different recordings having differ- 
ent mixtures of features and different HRV characteris- 
tics. For each test record and for each selected portion of 
the test record, the HRV parameters obtained from the 
commercial device should be compared with the known 
characteristics of the initial signal. .Any discrepancies 
found should be itemised in respect of features intro- 
duced Into the recording, e.g. errors caused by increased 
noise, errors due to fiducial point wander, etc. System- 
atic bias introduced by the equipment as well as its 
relative errors should be established. 
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Reporting the results 

A technical report of the testing should be prepared 
solely at the testing site independent of the manufac- 
turer(s) of the tested device. 
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